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Abstract 

Simple proofs of the Hermite-Biehler and Routh-Hurwitz theorems are presented. The total non¬ 
negativity of the Hurwitz matrix of a stable real polynomial follows as an immediate corollary. 


1 Introduction 

The classical result of Routh-Hurwitz on the stability of polynomials is now more than a century old. Its 
rich connections with other areas of analysis and algebra have been exposed in many subsequent works. 
Monographs by Postnikov m and by Chebotarev and Meiman (2 give a detailed account of such related 
questions, including the amplitude-phase interpretation of stability, Sturm chains, Cauchy indices, the prin¬ 
ciple of the argument, continued fractions, Hermite-Biehler theorem, and rational lossless functions. The 
interested reader should also see Barnett and Siljak’s centennial survey JJJ and references therein to find out 
what control theory problems can be solved using Routh’s algorithm. Krein and Naimark’s survey m and 
Householder’s article El explore connections of the Routh-Hurwitz scheme with Bezoutians, while Genin’s 
article |Sj emphasizes connections with Euclid’s algorithm and orthogonal polynomials, and presents a gen¬ 
eralized Routh-Hurwitz algorithm suitable, e.g., for testing nonnegativity of a polynomial. Asner (21 and 
Kemperman m found the link between stability and total positivity of the Hurwitz matrix. The Routh- 
Hurwitz algorithm, originally formulated for real polynomials, has been extended to complex polynomials 
(see [Zj or 0) and further to wider classes of analytic functions (see jjjj). This list is in no way exhaustive, 
since the existing literature on the subject is enormous. 

This note is not a survey of the field. Nor does it present a new approach to the subject. The note 
serves to derive, in a most elementary and economical way, three basic results in the Routh-Hurwitz theory, 
namely, the Hermite-Biehler theorem, the Routh-Hurwitz criterion, and the total positivity of the Hurwitz 
matrix. Because of this last issue, the setup is restricted to real polynomials. However, the proof of the 
Hermite-Biehler theorem extends verbatim to the complex case by considering the (generalized) odd and 
even part of a polynomial. My approach is minimalistic. For example, orthogonality of polynomials, rational 
lossless functions and the like are not discussed when the only fact needed is interlacing of roots. The point 
is not that these connections are unimportant, but that they are not needed for a quick and direct derivation 
of the basics of the Routh-Hurwitz theory. 

*On leave from Department of Computer Sciences, University of Wisconsin, Madison, Wisconsin 53706 U.S.A. Supported 
by Alexander von Humboldt Foundation. 
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Two papers, [I and m, offer alternative elementary proofs of the Routh-Hurwitz scheme (labeled 
Theorem |2 in this note). The proof in pQ is based on geometric considerations in the complex plane and the 
proof in d, simpler in my opinion, solely on continuity of the roots of a polynomial. Here, on the other hand, 
three results are derived essentially at once, Routh-Hurwitz being a direct consequence of root interlacing 
obtained in Hermite-Biehler, and the total nonnegativity of the Hurwitz matrix a direct consequence of 
Routh-Hurwitz interpreted as a matrix factorization formula. 


2 Proofs 


Definition. A polynomial / is stable if the condition /(z) = 0 implies Re 2 < 0. 

The following is a version of the Hermite-Biehler theorem (0i El)- 

Theorem 1 Let /( x) = p(x 2 ) + xq{x 2 ), f{x), p(x), q(x) £ IR[x]. The following are equivalent. 

A. The polynomial f is stable. 

/ 2 \ 

B. The polynomials p(—x 2 ) and xq(—x 2 ) have simple real interlacing roots and Re i) > 0 for some zq 

with Re zq > 0. 


Proof. A=^B: If / is stable, then /(z) = a J~[ ; (z — Zj ) with all Zj in the left half-plane. If Imz > 0, 
then |i z + Zj\ > | iz — Zj\ for all j , hence |/(iz)| > |/(iz)| or, by expanding the squares of both absolute values 
and simplifying, Imp(—z 2 )zg(—z 2 ) < 0. This implies that the functions 


p(-z 2 ) zq(—z 2 ) 

Z t—> -, 0 - 

zq(—z 2 ) p{—z 2 ) 

take on real values only on the real axis. Hence any non-trivial real linear combination 


Cl) 


\p(-z 2 ) + pzq(-z 2 ), A 2 + /i 2 ^0, 


( 2 ) 


has only real roots. Next, gcd(p(—x 2 ), xq{—x 2 )) = 1, for if not, then / would have either two roots with 
opposite real parts or one on the imaginary axis. So, if 0 had a multiple root, one of the functions 0 
would have a high-order crossing with some horizontal line. But if g(x ) — r = (x — xq) k h(x), h{x o) ^ 0, 
for analytic functions g , /i, and k > 1, then, for small e > 0, the equation g{x) = r — e k has solutions 
Xq + e 17r ^ 1+2 - , l/ /c /i(xo) _1 ^ fc e + o(s), j = 1,..., k. Hence the function g takes on real values somewhere off the 
real axis. This shows that no combination 0 has a multiple root. This also implies that the roots ofp(— x 2 ) 
and xq{—x 2 ) interlace, for if not, then one of the functions would preserve its sign on the interval between 
two consecutive roots of the other, hence, by a standard argument, there would be a combination 0 with 
a multiple root inside that interval. 

B=>A: If Condition B holds, then the function z i—> Re ^ zq(^) ) d° es n °t change its sign in the half-plane 

Rez > 0 and that sign is positive. So, the equation +1 = 0 or, equivalently, /(z) = 0, has no solution 
with Rez > 0. The roots of p(—x 2 ) and xq(—x 2 ) are distinct, so there is no solution to /(z) = 0 on the 
imaginary axis either. □ 


Remark. The beginning of this proof is in the spirit of the argument from pp. 13-15]. demonstrating 
that A implies that 

/ p (- z 2 ) \ 

Im ( — -— Im z < 0 whenever z ^ 1R. 

' 9 (~)/ 


The following Theorem is the essence of the Routh-Hurwitz scheme. It is proved in monographs using 
Cauchy indices, Sturm chains or the principle of the argument (see, e.g., 0 pp. 225-230]). A nice elementary 
proof is given in ni Here is a different elementary argument based on Theorem 0 
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Theorem 2 The polynomial f(x) = p(x 2 ) + xq(x 2 ), with p{x), q{x) G ]R[x], is stable if and only if 
c:=p(0)/g(0) > 0 and the polynomial /( x) = p{x 2 ) + xq(x 2 ) is stable, where p{x) :=q(x), q(x) := ^(p(x) — 
cq{x)). 

Proof. Necessity. Condition B in Theorem Q] is equivalent to p and q satisfying p(0)q(0) > 0 and 
having only simple zeros, all negative, interlacing, the rightmost zero being that of p. 

Let the pair (p, q) satisfy Condition B, let x n < ■ ■ ■ < x\ be the zeros of p , and yk < ■ ■ ■ < Ui the zeros 
of q , and assume wlog that p(0) > 0. Then, with y n any point to the left of x n in case k = n — 1, we have 
p and q of opposite sign in (yj .. Xj ), all j, and also (—1 ) 3 p(yj) > 0 for all j. But then, for any c > 0, the 
polynomial r := p — cq has the same sign as p on [yj .. Xj ), all j. In particular, also (—1 ) 3 r(yj) > 0, all 
j, and this implies that, in each of the n — 1 intervals (jjj+i --yj), r has an odd zero. If now, specifically, 
c = p(0)/g(0) (which is positive, by assumption), then r also has a zero at 0, and since its degree is no 
bigger than n, those n — 1 odd zeros must all be simple. But this implies that q is of degree n — 1, with 
all its zeros simple and negative, and these zeros separate those of p := q , and, in particular, q(yi) has the 
sign opposite to r(y±), i.e., to p(y i), i.e., is positive, hence both p and q are positive at 0. In short, if (p, q) 
satisfies Condition B of Theorem 0 then so does the pair (p, q). 

Sufficiency. Suppose f(x) is stable and c > 0. Since p(x 2 ) = cp(x 2 ) + x 2 q(x 2 ), q(x 2 ) = p(x 2 ), and, by 
Theorem ^ Re ( j- ) > 0 whenever Re^ > 0, one obtains 

\ zq(z 2 )J 

Re('h^_l = Re(£)+Re(^)>Ite0>O whenever Ree>0. 

Finally, if p[x 2 ) and xq(x 2 ) are relatively prime, so are p{x 2 ) and xq(x 2 ). This proves that Condition B of 
Theorem Q is met. □ 


Theorem [5] implies the following version of the Routh-Hurwitz theorem (jTB:, jn] ). 


Theorem 3 The polynomial f(x) =: ^”_ 0 ajX J with ao > 0 is stable if and only if its infinite Hurwitz matrix 
H{f) is a product of the form 

H(f) = J(d) ■ • • J(c n )H(b), (3) 

with all parameters Cj, j = 1,..., n, positive, and b a positive polynomial of degree 0. Here 
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Proof. Prove by induction that a polynomial / of degree n, /(0) > 0, is stable if and only if the 
first n + 1 leading principal minors A j(f), j = 1,... ,n + 1, of H{f) are positive and the factorization @ 
holds. Indeed, if n = 0, both properties are valid trivially. If n > 0, then, by the Lemma, / is stable if 
and only if c > 0 and / is stable. But, as one readily verifies, H(f) = J(c)H(f), hence, in particular, 
Aj+i(/) = c /(0)Aj(/), j = 0,1,...; here Ao:=l. Since deg/ = deg/ — 1, / satisfies the inductive 
hypothesis, hence so does /. □ 


Theorem 4 The Hurwitz matrix of a stable polynomial f satisfying /(0) > 0 is totally nonnegative. 

Proof. By Theorem [3 the factorization m holds with all parameters positive. By inspection, each 
factor is totally nonnegative, hence their product H{f) is also totally nonnegative. □ 

Theorem 2] was first proved in [5] and m 
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